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The effect of electron beam pitch angle and density gradient on solar type III radio 
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1.5D Particle-In-Cell simulations of a hot, low density electron beam injected into magnetized, 
maxwellian plasma were used to further explore the alternative non-gyrotropic beam driven elec- 
tromagnetic emission mechanism, first studied in Ref. PJ. Variation of beam injection angle and 
background density gradient showed that the emission process is caused by the perpendicular com- 
ponent of the beam injection current, whereas the parallel component only produces Langmuir 
waves, which play no role in the generation of EM waves in our mechanism. Particular emphasis 
was put on the case, where the beam is injected perpendicularly to the background magnetic field, 
as this turned off any electrostatic wave generation along the field and left a purely electromagnetic 
signal in the perpendicular components. The simulations establish the following key findings: i) 



Initially waves at a few tOce/j are excited, mode converted and emitted at 



ii) The emission 



intensity along the beam axis is proportional to the respective component of the kinetic energy of 
the beam; iii) The frequency of the escaping EM emission is independent of the injection angle; iv) 
A stronger background density gradient causes earlier emission; v) The beam electron distribution 
function in phase space shows harmonic oscillation in the perpendicular components at the relativis- 
tic gyrofrequency; vi) The requirement for cyclotron maser emission, -g£- > 0, is fulfilled; vii) The 
degree of linear polarization of the emission is strongly dependent on the beam injection angle; viii) 
The generated electromagnetic emission is left-hand elliptically polarized as the pitch angle tends 
to 90°; ix) The generated electromagnetic energy is of the order of 0.1% of the initial beam kinetic 
energy. 

PACS numbers: 



I. INTRODUCTION 



Solar type III radio bursts are believed to be a result 
of relativistic electrons of ~ 10 — lOOkeV, that propagate 
through solar plasma. The electrons might be accelerated 
in impulsive solar flares as they are observed to occur in 
groups of 10 and more 0. Potential sources for accelera- 
tion include magnetic reconnection and dispersive Alfven 
waves [H, HI i which are both likely to be found in flare re- 
gions. 

Observationally, the electromagnetic emission spectrum 
of type III bursts show a characteristic drop of frequency 
over time in the 'dynamical spectrum' (frequency versus 
time intensity plot). The frequency of the observed elec- 
tromagnetic signal typically drops a few orders of mag- 
nitude within minutes. The emission occurs at two sepa- 
rate frequencies, the fundamental plasma frequency and 
its second harmonic. It is observed that type III bursts 
are accompanied by electrostatic Langmuir waves as well 
as evidence for electron distribution functions that show 
a bump in the forward direction (parallel to the back- 
ground magnetic field) in phase space (5|. Further, there 
is evidence for back propagating Langmuir waves as well 
as ion acoustic waves. 

The common theory suggests that fast electrons, that 
propagate outward on an open magnetic field line 
through the Sun's atmosphere, trigger the 'bump-on-tail' 
instability, which is responsible for the generation of the 
observed Langmuir waves. The Langmuir waves then, 
in turn, generate electromagnetic emission at the local 



plasma frequency and its second harmonic. As the beam 
of fast electrons propagates from denser into less dense 
regions, the frequency of the generated emission drops, 
because the plasma frequency is a function of the square 
root of the density. There are various theories on how the 
electromagnetic emission is generated. Those theories in- 
clude non-linear wave-wave interactions (interactions of 
Langmuir waves with backward propagating Langmuir 
waves or ion-acoustic waves) [g, [7[, linear mode conver- 
sion on density gradients (partial reflection and partial 
conversion of a Langmuir z-mode wave)Q, and the an- 
tenna mechanism (trapping of Langmuir waves in de nsity 
cavities to drive currents at the second harmonic) [9l [Tof. 
Recent works [j], [TTJ] studied type III burst scenarios in 
particle-in-cell (PIC) type simulations. Ref. [ll| inves- 
tigated a magnetic field setup near a magnetic recon- 
nection region in a 2.5D (2D in space, 3D in velocities, 
fields), Ref. (l| considered a generic magnetic field line 
connecting Sun to earth in a 1.5D. Simulations in Ref. 1] 
were carried out with EPOCH, a fully relativistic, elec- 
tromagnetic PIC code. In the present study of solar type 
III radio bursts, EPOCH was also used to simulate dif- 
ferent beam injection angles into magnetised plasmas, as 
well as the effect of different background density gradi- 
ents on the generated emission. The simulations were 
set up to inject a mildly relativistic, hot, low density 
beam into a spatially ID plasma, that featured a con- 
stant background magnetic field as well as a parabolic 
density profile, both parallel to the spatial axis. The de- 
fault setup and analysis tools were tested by successfully 
reproducing the results published in Hef. [jj. Ref.[l| es- 
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tablished that generation of electromagnetic emission via 
a non-gyrotropic beam injection is possible. The primary 
goal of this study is the effect of different beam pitch an- 
gles and density gradients on the electromagnetic wave 
generation, as well as proving that EM emission can be 
generated without invoking Langmuir waves or the classi- 
cal plasma emission mechanism, which involves Langmuir 
and ion-sound waves satisfying resonant beat conditions 
to mode convert into escaping transverse EM radiation. 
In order to investigate the effect of different beam injec- 
tion angles, components of the beam injection momenta 
were varied, while the total beam momentum (and hence 
the total kinetic energy of the beam) was kept constant 
in all cases. Results showed that a parallel beam momen- 
tum resulted in generation of an electrostatic Langmuir 
wave, visible in the parallel electric field component, trav- 
elling at the speed of the beam (essentially Langmuir tur- 
bulence following the beam) . There was no excitation of 
Langmuir waves found, when the beam is injected per- 
pendicularly to the background magnetic field. In the 
perpendicular components of electric and magnetic fields, 
waves, that travel approximately with the speed of light, 
were found. This electromagnetic signal was also found 
for a perpendicular beam injection. Hence, the simulated 
generation of electromagnetic emission seems indepen- 
dent of the generation of electrostatic Langmuir waves. 
It should be noted, that, in 1.5D geometry, the classic 
plasma emission mechanism is switched off, because it 
involves interaction of Langmuir, ion-acoustic and elec- 
tromagnetic waves, that, when treated mathematically, 
involves a factor of |1cl X kx| 2 , with ki, being wave vector 
of the Langmuir wave and kp for the transverse electro- 
magnetic wave. The angle between the two wave vectors 
for the 1.5D case is equal to zero, setting also the inter- 
action probability to zero. 

In the case of perpendicular beam injection, the electron 
distribution function will not show a bump in the forward 
direction, therefore, the instability, that is supposed to be 
responsible for the Langmuir wave generation, is never 
created. This result suggested an alternative emission 
mechanism. 

Closer analysis of the distribution function in (parallel 
and perpendicular) momentum space showed, that, for 
non-parallel beam injection, a positive gradient was to 
be found in the perpendicular direction. Such shapes of 
distribution functions are referred to as cyclotron maser 
or loss-cone distributions. In special cases, if distribution 
functions are shaped appropriately, they are also called 
shell or horseshoe distributions, as studied in Ref . [l2l. Il3j 
in the context of auroral kilometric radiation (AKR). It 
is known that these distributions allow for wave genera- 
tion in most prominently x- and o-mode, but also z-mode 
[T3|. Which mode is growing strongest is prescribed by 
the ratio w ce /w pe , i.e. how the magnetic field compares to 
plasma density, where oj ce is the electron gyrofrequency 
and uj pe is the plasma frequency. 

This paper is structured as follows: In Sec|TT]we will dis- 
cuss in detail the simulation setup. We will vary the 



beam injection angle in Sec lIIII and demonstrate the im- 
pact of different background density gradients in SeclIVI 
In Sec[V] we will study the behaviour of the distribution 
function, while, in Sec lVIl we will analyse the polariza- 
tion of emitted waves as well as relate the generated field 
energy to the initial beam kinetic energy. Finally, we will 
summarize key points in IVIIl 



II. SIMULATION SETUP 

All simulations, presented in this paper, use EPOCH, 
a fully electromagnetic, relativistic particle-in-cell code 
that was developed by the Engineering and Physical Sci- 
ences Research Council (EPSRC)-funded collaborative 
computational plasma physics (CCPP) consortium of UK 
researchers. 

We are considering a 1.5D maxwellian plasma, allow- 
ing spatial variation in x only while electromagnetic 
fields and particle momenta have all three components. 
The background magnetic field along this line is kept 
constant B = B x = 0.0003T= 3G, setting the elec- 
tron gyrofrequency to Lo ce — 5.28 x 10 7 Hz rad every- 
where. The background temperature is T = 3 x 10 5 K 
and isotropic. The background plasma density at x = 



IS Hq 



10 14 m 3 , giving 



ui 



5.64 x 10 8 Hz rad. This 



sets ^ = 0.0935 < 1. Further, the electron Debye 
length is \ De = 3.78 x 10 _3 m. 

The simulation setup is such that we simulate a sin- 
gle magnetic field line connecting Sun and earth, while 
the grid size is Xjj e /2. It is predicted that in the limit 
r >• Rq the plasma density n e (r) cx r~ 2 ([THEBl), there- 
fore, we make use of 



n e .i(x) = n 



~\~ Tl-\- 



(i) 



with x m ax being the total system length and 



n + 



io- 



(2) 



resembling a parabolic density drop off, such that 
n e .i(0) = n e .i(x max ) = n and n e<i {x max /2) = n_n . 
The unrealistic growing density in the region x max /2 < 
x < x max is employed in order to use periodic boundary 
conditions. When analysing results we, therefore, draw 
conclusions from the region < x < x max /2 only. 
Further, in section ITVl we consider modified plasma den- 
sity profiles to study the effects of different density gra- 
dients. The following density profiles are used to study 
a weak and strong gradient 
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density profiles 



FIG. 1: The electron beam (solid), weak (dotted), moderate 
(dashed), strong (dotted-dashed) density gradients at t = 0. 
Densities are normalized to their maxima. 



n w {x) = n 
n s {x) = n 



•Kmax /2 ~\~ 

X Xmax /2 
Xmax I ^2 "T" 



n_ = 0.5 (3) 
n_ = 10~ 8 .(4) 



A beam of accelerated electrons is injected at sim- 
ulation time t — 0. It carries a total momentum of 
Pb = m e 7| with 7 w 1.155. The beam temperature is 
T b = 6 x 10 6 K. The peak beam density is n b0 = 10 n m~ 3 , 
while its spatial profile is defined as 



rib{x) 



riboe 



-[(x-x max /25)/(x ma!c /40)f 



(5) 



setting the beam density maximum at x max /25. In 
order to demonstrate the situation at t = in terms of 
densities as a function of space, we plot the beam and 
background number density profiles in FigJTJ Note that 
the beam is injected only once, at t = 0, i.e. the beam 
electrons are not replenished. 

As a result of the above defined quantities, at x = 0, 
the plasma beta is j3 — 0.0115. The mass ratio used is 
m,i/m e — 1836. 

Each simulation is run on 32 Dual Quad-core Xeon pro- 
cessors with 4 Tb of RAM. The simulation makes use of 
65000 grid points and 10000 particles per cell in the case 
of the background electrons and ions. Despite the low 
beam electron density (n ey i(0)/nb = 10 3 ), the simulation 
achieves a reasonable number of beam electrons per cell, 
since the beam is focussed on only roughly 10 simulation 
cells, resulting in 100 beam electrons per cell. Herein lies, 
however, a shortcoming of the simulation setup as real- 
istic background-to-beam electron ratios of « 10 5 — 10 7 
cannot be achieved, because they would require consid- 
erably larger computational resources. 



III. VARYING BEAM INJECTION ANGLE 

In this section we analyse the effect of a variation of 
the beam injection angle (with respect to the background 
magnetic field) on the EM emission. We change the in- 
jection pitch angles via variation of a perpendicular com- 
ponent in the beam momentum, p by ~ pf,sin#, while 
keeping the total momentum p b constant in all cases. 
Therefore, an increase of the beam pitch angle 9 is equiv- 
alent to an increase of perpendicular beam momentum 
Pb y , respectively a decrease in parallel beam momentum 
Pbx — Pb cos 8, hence, a decrease in propagation speed. 
For the case 9 = 90°, the parallel component of the beam 
momentum vanishes and the beam is trapped at its point 
of injection. The background plasma density profile in 
this section is always represented by Eq.([T]) which cor- 
responds to the dashed line in Fig|T] We study cases of 
9 = 15°, 45°, 60°, 90°. In order to analyse the simula- 
tion data, we make use of time-distance plots and spatial 
wavelet transforms of the electromagnetic field compo- 
nents. Time-distance plots offer the possibility of study- 
ing the development of a physical quantity with respect 
to space and time. We make use of this technique to 
plot components of electric and magnetic fields as well 
as the electron density. We observe waves in both elec- 
tric and magnetic components. The slope that is be- 
ing produced by the wave front in this plot can then 
shed light on the propagation speed of that wave. We 
choose units for spatial and temporal axis such that a 
slope of « 1 corresponds to propagation with the speed of 
light. Hence, finding waves that travel with w c suggests 
that the emission is electromagnetic. Dynamical spec- 
tra show a characteristic feature of solar type III bursts. 
These plots show the development of the emission fre- 
quency over time. In the case of solar type III bursts, the 
dynamical spectrum shows a frequency drop over time. 
This drop is due to the decreasing background plasma 
density, that a beam of accelerated electrons encounters, 
while travelling away from the Sun. The drop in den- 
sity correlates to a drop of the local plasma frequency, 
since ui pe oc ^fn^. As the beam triggers emission at the 
local plasma frequency (and its second harmonic), while 
it travels from denser into less dense plasma, it generates 
emission of dropping frequency that is characteristic for 
solar type III bursts. We investigate the frequency of the 
generated emission by performing a wavelet analysis of a 
perpendicular component of the electric field; we choose 
E y . We take a snapshot of E y at a given time and per- 
form a spatial wavelet analysis. The wavelet software was 
provided by C. Torrence and G. Compo[22j]. 

Figs l2l3l4l5l show time-distance plots of components of 
the electric field in the left column. The right column 
contains plots of changes in number density, the magnetic 
y-component and the wavelet transform of E y taken at 
t = lOOcjpg 1 . Electric and magnetic field strengths are 
normalized to units of uj pe cm e /e and tu pe m e /e, respec- 
tively, while distance and time are measured in c/uj pe 
and w" 1 , accordingly, where uj pe is the plasma frequency 
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FIG. 2: 6 — 15°. The moderate density profile, given by Eq.JT]), was used. Left column: time-distance-plots of electric field 
components. Right: time-distance-plot of changes in density, magnetic field y-component, and the wavelet transform of E v . 
Note that the background for the wavelet transform was set to white colour and does not refer to maximum amount of emission 
on the sides of the plot. Further, the black lines track the local plasma frequency and its second harmonic. 



at x = 0. This is a good value of reference, as the back- 
ground density (and therefore the plasma frequency) at 
this point are equal in all simulations, regardless of which 
density profile for the background (see FigQJ is used. 
Comparing the time-distance-plots for electromagnetic 
fields, we find wave generation in the perpendicular com- 
ponents for all injection angles. The wave front shows a 
slope of pa 1, is therefore travelling with the speed of light 
and can be interpreted as electromagnetic emission. The 
wave propagation speed in the perpendicular components 
does not depend on the beam injection angle and is there- 
fore independent of the beam propagation speed. Similar 
analysis of the parallel component E x shows rather differ- 
ent behaviour. In the 8 = 90° case (FigJS|), E x does not 
show any wave features and no parallel component to the 
emission is generated. In the other cases (Fig l2l3l4")l . we 
can see generation of a standing wave as well as a travel- 
ling electrostatic wave, which shows different propagation 
speeds for each of the cases, and can be associated with 
the beam propagation speed. This signal also has com- 



ponents in the perpendicular plane, most notable in the 
magnetic B y component, but also in (rather weakly) E y 
and E z , where it interferes with contributions of the elec- 
tromagnetic emission. The width of the signal at t = 
corresponds to the width of the beam at the start of the 
simulation. 

In all cases, the density profile remains unchanged (see 
top right panels on Figsj2][5]). This is important as it ex- 
cludes any possibility of developing density cavities in the 
plasma, which could trap Langmuir waves and contribute 
to the emission in form of the antenna mechanism [loj . 
Thus, in our simulations, EM emission is generated by a 
mechanism other than the antenna mechanism or plasma 
emission (the latter requires at least 2D in space). 
Comparing wavelet transforms for various pitch angles, 
we clearly see that emission intensity is proportional to 
the perpendicular component of the beam momentum. 
Other than that, there is no evidence of change in the 
wavelet transform and emission characteristics as peak 
frequencies, emission times (respectively location) and 
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wave packet shapes are all preserved. Analysis of the 
peak emission intensity as a function of pitch angle is 
shown in FigEl 

The figure shows that the peak emission intensity (nor- 
malised to its maximum value at = 90°) can be al- 
most perfectly approximated by sin 2 (9), and therefore 
pl y . This result suggests a direct correlation of the beam 
kinetic energy, E\ in oc p\ and the generated emission 
intensity. 

IV. TIME DEVELOPMENT AND VARYING 
THE BACKGROUND DENSITY PROFILE 

We now focus on the 8 = 90° case, as it shows the 
maximum generated emission in E y . In this section we 
investigate the role of the background density profile. We 
refer the reader to FigQ]for the different density profiles. 
Additionally, we investigate a case without any density 
drop, where n c [x) = uq. 

The time development of the wavelet transform allows 
insight into the emission generation process, wave prop- 
agation, as well as frequency drift. We generated movies 



to demonstrate this (movie 1, movie 2, movie 3, and 
movie 4 in Ref.[I3|). Again we focus on the wavelet 
transform of E y . We include curves that mark the lo- 
cal plasma frequency and its second harmonic, as well as 

uj z = — |w ce + ^^Juj^ e + 4ojp e , which is the cut-off for the 

z-mode. From the movies, it is evident that as soon as 
the beam is injected, we observe a pulsating signal of gen- 
erated emission. The frequency of this emission is much 
smaller than oj pe and, therefore, does not allow instant 
propagation (note that uj ce /ujp e <C 1). Instead, pulsation 
eventually dies down only before an apparently stable 
wave packet is formed. The stable packet then drifts 
to higher frequencies. This suggests mode coupling of 
the cyclotron emission on the gradient to (possibly) a 
z-mode. When the peak frequency of the wave packet 
is comparable to the local plasma frequency (which for 
our parameters is not significantly different from w z ), the 
packet starts propagating. Interestingly, the increase in 
frequency does not stop at this point but only when the 
wave packet travels into a region where the second har- 
monic of the local frequency, 2fi oca i, matches the fun- 
damental of the plasma frequency at the injection point 
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FIG. 4: As FigE] but 6 = 60°. 



finj- (Note that the beam is not injected at x = 0, but 
has a finite width and fi n j is the plasma frequency at the 
beam density function's maximum at x maa: /25.) From 
this point on the wave packet moves on with constant 
frequency. 

In FiglTl we show snapshots of wavelet transforms at dif- 
ferent simulation times. We choose to look at an early 
snapshot of the pulsation phase, a second snapshot when 
the wave packet couples to frequencies « f pe and, ulti- 
mately, the final snapshot in the simulation. Note that 
for the mid row case, the total simulation time has been 
extended to t = 150o;~ e . In order to not get any inter- 
ference due to periodic boundary conditions, the total 
system size was increased by a factor of 1.5 in this run. 

We see that for constant density, there is no emission, 
as the wave packet never mode couples to frequencies that 
would allow propagation. At this point, we sum up the 
key findings: The beam injection evokes pulsating emis- 
sion far below the local plasma frequency f pe . Eventu- 
ally, this cyclotron emission mode couples to the z-mode, 
which has a frequency comparable to the local plasma fre- 
quency, and that allows for escape of the emission. The 
wave packet frequency increases up to the point where 



it reaches the plasma frequency of the original beam in- 
jection point, fi n j. This happens at the point where the 
second harmonic of the local plasma frequency, 2// oca /, is 
of the order of fi n j . 

These observations enable us to draw the following con- 
clusions with regards to the role of the density profile in 
our simulations: In case of the constant density profile, 
there is no escaping emission, as the generated emission 
cannot couple to a propagating wave mode. A steeper 
gradient in the density profile will result in earlier emis- 
sion, because the mode coupling process couples to a fre- 
quency that is able to escape sooner than it does for 
weaker gradients. 

We would like to note that Fig[7] as well as all movies 
of wavelet transforms in Ref . [17| suggest that the pulsa- 
tion maximum is not at the y-axis value corresponding 
to f ce . This is due to the fact that in order to be able to 
interpret the y-axis as frequency, one needs to be able to 
relate space and time via x — ct. The latter is only appli- 
cable to escaping electromagnetic radiation. In the left 
bottom corner of the wavelet transforms, i.e. the beam 
injection region, this relation does not hold, as we are 
dealing with a trapped mode far below f pe . Therefore, 



7 




FIG. 5: As FigE] but 9 = 90°. 



whilst the wave packet resides in this region, the y-axis 
actually corresponds to the wave number k rather than 
the wave frequency. Numerical runs of wider (respec- 
tively narrower) beam width, not shown here, corrobo- 
rate that the y-value of the spatial wavelet transform in 
regions well below f pe scales as A:. It is only after / w f pe 
and x — ct applies, that the y-axis can be interpretted as 
wave frequency. Further, the fact that the initial pulsa- 
tion is excited at the relativistic cyclotron frequency can 
be evidenced from figures equivalent to Figs l2l3l4l5l but 
for longer runs, which are not shown here. E y pannels 
of the figures shown here already suggest that there is a 
time interval of w 60WZ, 1 between the bottom left wave 
structure and the one following. Thus, the frequency of 
the oscillation is 1/At = 27r/ pe /60 « 0.1/ pe w f ce . 

V. DISTRIBUTION FUNCTION DYNAMICS 

A. Cartesian Coordinates 

At any given timestep, in EPOCH, distribution func- 
tions for every species are functions of space and momen- 



tum, i.e. f(x,p). They satisfy the equation 

J f (T {x,p)d 3 pdx = N a (6) 

with N a being the total number of particles of a species 
a (electrons, ions, beam electrons). The distribution 
functions of both the background plasma and the beam 
electrons in momentum space are initially maxwellian. 
This allows us to separate the spatial from the momen- 
tum part, by f(x,p) = f(x)f(p), at t = 0. The momen- 
tum part can, then, be described by 

fe.Mp) = h^e-^y+^y^^ (7) 
for the background and 

/ 60 (p) = n be -l(Px-Pbx) 2 + (Py-Pby) 2 +pl]/<,2mek B T b ) fg\ 

for the beam respectively, where ks is Boltzmann's 
constant and n e ,i,& are normalization constants corre- 
sponding to a species (electrons, ions, beam electrons). 
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FIG. 6: The maximum emission intensity lmax(0) normalized 
to Imax (90°) (crosses) and sin 2 ((9) (dashed) as a function of 
pitch angle 9 in degrees. 



EPOCH allows for distribution functions to be com- 
puted with respect to spatial and momentum coordi- 
nate components, as well as all desired combinations. It 
is therefore possible to generate arrays for distribution 
functions in all momentum components, f(j> x ,p y ,p z ), 
but also specific components, for example f{jp x ),f(p x )... 
Here, the spatial dependency has already been integrated 
out by EPOCH. Following the time development of the 
distribution function for various angles, we refer to movie 
5 and movie 6 in Ref . [I?} • In the movies we show the de- 
velopment of the distribution function in three panels, 
top to bottom: f{p x ), f(p y ), f(p z )- 

Fig J5] shows snapshots of f(p y ) and f(p z ) at different 
times. Each panel, generally, shows two bumps. The 
larger one refers to the background electrons, while the 
smaller one corresponds to the beam. We can clearly see 
that the background distribution function as well as the 
beam distribution in parallel momentum space (B along 
x) remain almost constant throughout our simulations. 
In perpendicular momentum space, however, we observe 
oscillation of the beam distribution, as well as a broad- 
ening and widening of the beam bump. We focus on the 
perpendicular panels and follow the position of the beam 
maximum in time, which is shown in Fig|9l 

FigE] is representative for all cases studied. The max- 
ima for this run are at % « 0.5775r7i e c. (Note that 
for different injection angles, the maxima are shifted to 
different values on the ordinate, because different injec- 
tion angles result in different values of p y , but the func- 
tional characteristics will stay the same.) It shows that 
the maximum of the beam distribution function oscil- 
lates in the perpendicular components. The frequency 
of the oscillation is the relativistic cyclotron frequency, 
uj cr = (Jce/l- It is trivial to deduce that the beam is gy- 
rating around the magnetic field. The evident gyration 
along with the observation, that the initially generated 
pulsating emission is of the order of the cyclotron fre- 



B. Parallel and Perpendicular Coordinates 

While the cartesian form offers a valuable look into 
dynamics of the beam, investigation of instabilities often 
require a transformation of the distribution function to a 
momentum space with coordinates parallel and perpen- 
dicular to the background magnetic field, i.e. f(p\\,px)- 
While in the more conventional plasma emission process, 
Langmuir waves are being generated via the 'bump-on- 
tail'-instability, referring to a bump in the distribution 
function in parallel momentum space, > 0, the re- 
quirement for the triggering of cyclotron maser (or loss- 
cone) emission is a positive slope of the distribution func- 
tion with respect to the perpendicular direction, > 0. 
EPOCH does not allow for direct output of f(p\\,p±), so 
a transformation algorithm was developed. The trans- 
formation relations are straight forward 



P± 



pI+pI 



(9) 



The data volume accumulated by the simulation is 
prescribed by the resolution of the distribution func- 
tion. A resolution of 300 was chosen for an interval of 
— 0.91rn e c < p x , y , z < 0.91m e c, to retain reasonable reso- 
lution while keeping data volume comparatively small. 
The algorithm was tested by successfully reproducing 
Fig.l in Ref. [3]. 

Movie 7 in Ref. 17] shows the time development of the 
distribution function for the 9 = 90° case. The top left 
panel shows a plot of the distribution function on a log- 
arithmic scale at a given time step, while the top right 
panel plots the relative change with respect to the initial 
distribution, ^fejp~ (where e is a small number to avoid 
division by zero), in order to visualise dynamics. In the 
lower row, the gradients (left) and (right) are 
shown. We can clearly see that the initial distribution 
fulfils the criterion > 0. Further, the movie shows 
interesting dynamics in the shape of f(p\\,Px)' The left 
panel shows no drastic changes throughout the simula- 
tion time, while the relative change on the right panel 
suggest shifting to higher and lower perpendicular mo- 
menta. First, the distribution increases in momenta that 
lie below the initial maximum, 0.5775m e c, then the 
movement is reversed. Until around t w 30w pe we can 
see frequent change in the direction of shifting. Com- 
paring this to movie 3 in Ref.[l7|], we see that the time 
scales correspond to the pulsating phase in the wavelet 
transform, suggesting an energy exchange between beam 
and EM fields, as found in Ref.[l8j]. Eventually the dis- 
tribution shifts to momenta above the initial maximum 
and forms almost a straight line at w 0.5775m e c, sep- 
arating regions of increased distribution from decreased 



9 



t= 1 .6 4) " 1 t=50 i) " 1 t=100 u ~ 1 




* [co^ n ] k [my '] * [ciy n ] 



FIG. 7: Wavelet transforms for different density gradients at different times. Top row: density is constant. Middle row: 
moderate gradient according to dashed line in Fig[T] Bottom row: strong gradient according to dotted-dashed line in Fig[T] 
Black curves correspond to local plasma frequencies and the second harmonic, the horizontal line is the plasma frequency 
corresponding to the point of the beam density maximum at the time of the beam injection, finj- 
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FIG. 8: Distribution function for the background and beam 
electrons in momentum space at different times (from top to 
bottom) for — 90°; f(p y ) (left column) and /(p z ) (right 
column) 



FIG. 9: Position of beam maximum in momentum space as 
function of time for 6 — 90°; p y (solid) and p z (dotted) in 
units of [m e c]. 
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FIG. 10: Top left: f(p\\,pj_,i). top right: ^"{^ . bottom 

left: MS. bottom right: at t = 40^. Gradients 

are shown on logarithmic scales, which cannot show negative 
values, therefore, they appear white as the background. 



with the initial bulk momentum < po >, given by 
fpf (x,p)d 3 p 



< Po > 



n(x) 

fo(x) J pf Q (p)d 3 p 
n(x) 



Using the relation 



/ 



f(p)d 3 p = 1 



(13) 



f(x,p)d 3 p = n(x) (14) 
along with the normalization condition 



(15) 



we see that the spatial part of the distribution function 
becomes the particle density distribution 



ones, at t pa 40w pe . Comparison with movie 3 in Ref.[17j| 
shows that this is about the time, when we can see a sta- 
ble (no longer pulsating) wave packet being formed. At 
t pa (70 — 90)u>pe, the distribution function shifts again 
to momenta below pa 0.5775m e c. The wavelet transform 
shows that this is the time, when the wave packet passes 
over the local plasma frequency limit. Fig[TU]is a repre- 
sentative snapshot at t — 40w pe 1 . The gradients in the 
bottom row are presented on a logarithmic scale, there- 
fore, negative gradients cannot be shown and appear in 
the same white background colour. 

VI. GENERATED ELECTROMAGNETIC FIELD 
ENERGY AND POLARIZATION 

We calculate the electromagnetic field energy by use of 
the Poynting theorem, 

w{t) = f [^e E(x, tf + ^-B(x, t) 2 ]dx (10) 

We take care to exclude contributions from the back- 
ground magnetic field Bo. We relate the field energy to 
the initial kinetic energy of the beam, E^ m (0), which we 
calculate analytically using the expression for the mean 
rclativistic kinetic energy 

< E b k Z m (0) >= y"n b (a;)TO e c 2 [7(< p >) - l]dx (11) 

where 7(< po >) is the Lorentz factor in terms of 
initial bulk momentum 

7(< Po >) = Jl + (^^) 2 (12) 



fix) = n(x) (16) 
Hence, the bulk momentum can be written as 

< Po >= J pfo(p)d 3 P (17) 

where /o(p) is the initial beam distribution function 
and can be deduced from Eq.®. For our simulation 
parameters, the initial kinetic energy of the beam is 
E b k e ^ m (0) = 7.89272 x 10~ 3 J. We use this result to see 
how much of the beam energy gets converted into mag- 
netic energy, by plotting this in Fig llll 

The figure shows the time development of the ratio 
ui(£)/.E^ m (0) for different beam injection angles. We 
can see that the mechanism has a typical efficiency of 
pa 10~ 3 -10~ 2 . Most notable is that, for our parame- 
ters, a smaller injection angle 9 results in a higher ef- 
ficiency. Note that this is not in contradiction with 
Fig|51 Here, we take all electromagnetic field compo- 
nents into account, whereas, in Figj6] we focused solely 
on the peak intensity of the electric y-component, namely 
the wavelet transform of E y , and showed a direct re- 
lation to the y-component of the kinetic energy of the 
beam. Fig|6]is also in line with the previously presented 
time-distance plots Figsl^JSJ Comparing the scales in 
the respective color bars shows that the parallel compo- 
nent of the electric field, hence the electrostatic wake of 
the beam, is dominant for the presented cases (except 
of course for the 9 = 90°-case, where the parallel com- 
ponent vanishes). Fig llll shows the pulsating nature of 
the early stages of the emission generation. The pulsa- 
tion can be associated with the electromagnetic part of 
the emission, because it is generated in all cases. In the 
9 = 90°-case (FigJTTJi) , the initially generated electro- 
magnetic emission stabilizes at pa 10 _3 i?j^ m (0). In the 
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FIG. 11: Total field energy generated, w, normalized to initial 
kinetic energy of the beam, -Ejci^ m (0) according to Eq llll as 
a function of time for a) 9 = 15° b) 9 = 45° c) 9 = 60° d) 
6 = 90° 
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FIG. 12: Degree of linear polarization L as a function of time 
for a) 9 = 15° b) 6 = 45° c) = 60° d) 9 = 90° 



FIG. 13: Polarization in the perpendicular (y,z-)p\a,ne as a 
function of time for 9 — 90°. The :r-axis points out of the 
paper plane. The black 'x' marks the starting point, color 
corresponds to time in units of uj^} . Electric field components 
are given in units of Lj pe cm e /e. 

case. Not surprisingly, we see a very strong degree of 
linear polarization for this case. 

Concentrating on 8 = 90°, we plot the time evolution 
of the electric field in the perpendicular (?/,z-)plane in 
FigH3l In the figure, the a:-axis points out of the plane. 
The time development is then shown by a change of pixel 
color, therefore we see that the wave is left-hand polar- 
ized. 

A cyclotron maser in our parameter space (Lo ce /uj pe <C 
1) is expected to generate waves in the z-mode, at fre- 
quencies that are harmonics of the cyclotron frequency 
(l4| . The polarization of the z-mode is left-handed for 
uj < uj pe and right-handed for uj > oj pe (l9| . This is 
in line with our result, as the stable wave packet was 
generated below the plasma frequency, and is therefore 
expected to show left-hand polarization. 



other cases, we see that eventually the electrostatic wake 
of the beam overpowers the electromagnetic part signif- 
icantly and subsequently causes an increase in the total 
field energy up to » l(T 2 £;^ m (0) in the 6 = 15°-case 

(FigSBi). 

The degree of linear polarization L is defined similarly to 
Ref.pt 



W\\ + w± 

and its time development is presented in FigfT2l 
Fig[l2] shows a strong dependence of L with respect 
to the injection angle. As shown in FigfTTJ the paral- 
lel component of the field energy becomes dominant for 
small injection angles, but is negligible in the = 90°- 



VII. CONCLUSION 

In this paper, we study further the details of the 
non-gyrotropic beam driven emission mechanism, first 
outlined in Ref.[l[. We carried out 1.5D PIC simulations 
of a super-thermal electron beam being injected into a 
magnetised, maxwellian plasma. We investigated the 
role of the injection angle (pitch angle) of the beam with 
respect to the constant background magnetic field. 
The evolution of spatial wavelet transforms of E y was 
demonstrated in movies. The movies show that, initially, 
pulsating emission well below the local plasma frequency 
is generated. Waves of such low frequencies cannot 
propagate in plasma, thus, in order for the waves to 
escape, a mode conversion has to take place. In the 
movies, mode conversion is identified by the formation of 
a stable wave packet, which drifts to higher frequencies, 
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propagating only very slowly from the injection region. 
When the wave packet reaches frequencies of the order 
of the plasma frequency, propagation accelerates. The 
frequency drift, however, does not stop at this point, 
but only when the wave packet frequency reaches the 
frequency that corresponds to the plasma frequency at 
the beam injection point. Once this frequency is reached, 
the wave packet propagates without any frequency drift. 
Four cases of different density gradients were studied: a) 
a constant background density b) a weak gradient c) the 
default case as used by Ref.[l| and d) a strong gradient. 
The characteristics of the emission mechanism were 
not changed by a different density profile. However, a 
steeper gradient means that the wave packet can reach 
frequencies that allow for propagation earlier, while the 
weak gradient has the opposite effect. In case of the 
constant background, no emission could escape the beam 
injection region, as the necessary frequencies were not 
reached within the simulated time. Thus we conclude, 
that the mode coupling, namely the cyclotron emission 
coupling to the allegedly z-mode, is facilitated by the 
density gradient. 

We found that the intensity of the wavelet transform 
of E y was proportional to the y-component of the 
kinetic energy of the injected beam. There was no 
other significant influence of the variation of the beam 
injection angle 9. However, in the 9 = 90°-case, there 
was no electrostatic Langmuir wave signal found. This 
can be explained by the lack of a positive gradient 
in the distribution function with respect to parallel 
momentum. Nevertheless, generation of electromagnetic 
waves was evident. This sets the presented mechanism 
apart from conventional mechanisms, presented in the 
introduction. 

Further, movies were generated to show the evolution 
of the distribution function in phase space. On one 
hand, when plotted in cartesian coordinates, gyration 
of the beam around the magnetic field could be seen; 
on the other hand, when plotted on a plane of parallel 
and perpendicular momentum axes, positive gradients 
as well as subtle dynamics of the distribution function 
was observed. The key finding is that, indeed, the 
requirement for the cyclotron maser instability, > 0, 
is fulfilled throughout the simulation. Cyclotron masers 
in overdense plasmas (oj ce /u! pe -C 1) are expected to 
generate waves in the z-mode, at harmonics of the 
electron cyclotron frequency [lij • 

We would like to comment, why we think there is no 
continuous generation of cyclotron maser emission, 
despite the fact that > throughout the simulation. 
The quasi-linear relaxation, i.e. the plateau formation 
shown in Ref.[12|, that shuts off the instability, happens 
on the time scale of the inverse growth rate. The growth 
rate calculation of our parameter space will be presented 
elsewhere. The reason why the emission generation is 
not continuous is, that we do not replenish the beam, as 
it is only injected at t = 0. In solar flares, of course, the 
electron beam injection will have temporal (and spatial) 



extent. 

The distribution function shows subtle dynamics 
throughout the simulation. Some correlations can be 
seen between the evolution of the wavelet transforms and 
the distribution function, such as in the initial pulsation 
phase as well as at the wave packet ejection point. This 
needs to be investigated further. 

The degree of linear polarization of the escaping wave 
packet was strongly dependent on the injection angle. 
For the 9 = 90° left-handed polarization of the 

emitted wave was found. This is in line with z-mode 
wave generation, as the stable wave packet is actually 
formed below the plasma frequency and the polarization 
of the z-mode for u < uj pe is left-handed [l!| . At present 
we do not have definite proof that the generated emission 
is a z-mode. Such proof can be obtained by producing 
different numerical runs and producing (w, fc)-pairs and 
compare these to predictions of dispersion relations as 
in Ref.[19] which presents a comprehensive study, but 
not exactly for the same parameters, that are used here. 
The total generated electromagnetic field energy was of 
the order of 10 -3 - 10~ 2 of the initial kinetic energy of 
the beam and varied as the pitch angle was altered. 
The presented mechanism is relevant to solar type III 
radio bursts, as it is triggered by the perpendicular 
component of the electron beam, while the parallel 
component merely generates Langmuir waves. The 
present study focuses on the injection of a single beam 
as an initial value problem (i.e. without replenishing 
the beam), but in flares many beams are generated 
and due to the outward propagation, the cyclotron 
maser is triggered along the beam trajectory, which 
mode converts on the density gradient into the z-mode, 
thus generating electromagnetic emission of decreasing 
frequencies, that give the characteristic shape of the dy- 
namical spectrum for type III bursts. Ref.[20| indicates 
that beam pitch angles are a few degrees. This means 
that Langmuir waves would be generated, and yet, as 
shown in Ref . [2lj . Langmuir waves play no role in EM 
emission generation, because, in 1.5D, plasma emission 
is not possible. Reality is - of course - not 1.5D and, in 
higher dimensions, plasma emission will switch on; as 
well as the antenna mechanism, if density cavities are 
present. Additional studies are needed to measure rela- 
tive strength/importance of the different mechanisms in 
different physical situations. Further analysis may also 
help to better understand the mechanism presented in 
this study. Such analysis may include the determination 
of the mode of the escaping emission as well as detailed 
analysis of the mode coupling process. 

The authors are financially supported by the HEFCE- 
funded South East Physics Network (SEPNET). D.T.'s 
research is supported by The Leverhulme Trust Re- 
search Project Grant RPG-311 and STFC Grants 
ST/J001546/1 and ST/H008799/1. 

note added in proofs: After this work was complete, 
we became aware of the following: (i) when a case of 
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Tbackground = 3 x 10 8 is studied, eliminating the electron 
cyclotron maser (ECM) instability as < 0, Vp, EM 
emission identical to the one in this study still occurs; 
(ii) when a ring-shaped (in p y and p z ) beam distribution 
is considered, implying jj_ = 0, no EM emission is 
generated. These findings indicate that both effects - 
the ECM instability and EM emission from transverse 



currents - are present (and competing). It is probable 
that the ECM instability growth rate is so small that 
it cannot develop by the end of the simulation. Cal- 
culation of the growth rates commensurate to physical 
parameters of type III bursts and the ring distribution 
will be published elsewhere. 
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